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Abstract: Hawking radiation from black holes has been studied as a phenomenon 
of quantum tunneling of particles through their horizons. We have extended this 
approach to study the tunneling of Dirac particles from a large class of black holes 
which includes those with acceleration and rotation as well. We have calculated the 
tunneling probability of incoming and outgoing particles, and recovered the correct 
Hawking temperature by this method. 



1. Introduction 



Hawking's arguments [1, 2] for thermal radiation and evaporation of black holes were 
based on quantum field theoretic analysis in curved spacetime. This breakthrough 
gave rise to a whole new field of investigation where three different areas, the quan- 
tum theory, general relativity and thermodynamics, come together. A major devel- 
opment in this direction was to interpret Hawking radiation as a quantum tunneling 
phenomenon [3-5]. Using this procedure the Hawking temperature was calculated 
from the tunneling probabilities of incoming and outgoing particles [6, 7]. This gen- 
erated a lot of interest in studying radiations of different particles in the context of 
various black hole spacetimes. In this regard different semiclassical approaches were 
adopted to study radiation by scalar and Dirac particles. In particular, a procedure 
was set up and tunneling of charged and uncharged Dirac particles was studied for 
Kerr, Kerr- Newman, Taub-NUT and Godel spacetimes [8-10]. This approach was 
applied to other four-dimensional [11-15] and three-dimensional [16, 17] spacetimes 
as well. 

In this paper we study a large class [18-23] of Plebahski and Demiahski space- 
times for the tunneling procedure. We are basically motivated by two reasons. 
Firstly, the spacetime we study generalizes many black holes like the Schwarzschild, 
Kerr and Taub-NUT, and the tunneling probabilities and Hawking temperatures for 
these black holes can be recovered as special cases of our study. Secondly, the met- 
ric we study has very interesting interpretation in itself as accelerating and rotating 
black holes which have been studied in the literature for various properties. These 
black holes admit two rotation horizons and two acceleration horizons. In this pa- 
per we apply the procedure mentioned above to study tunneling of Dirac particles 
from these black holes. We calculate their tunneling probability and the Hawking 
temperature at these horizons. 



2. Accelerating and rotating black holes 

The metric for accelerating and rotating black holes without the cosmological con- 
stant can be written as [18, 19] 



,9 1 , ,Q a 2 P sin 2 9. , 9 p 2 , 9 p 2 , n9 

= — —^ + Q dr2 + -P de 

,P(r 2 + a 2 ) 2 sin 2 6 Qa 2 sin 4 ,,, 2a sin 2 0(P(r 2 + a 2 ) - Q)dtdA„ , , 
+ < ? — )W nW (2A) 
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where 

Q = l-arcos9, (2.2) 

p 2 = r 2 + a 2 cos 2 9, (2.3) 

P = l-2aM cos 9 + a 2 a 2 cos 2 9, (2.4) 

Q = (a 2 - 2Mr + r 2 )(l - a 2 r 2 ). (2.5) 

Here the parameters M, a and a represent the mass, angular momentum per unit 
mass, and acceleration of the source respectively. Following the notation of Ref. [9] 
the above metric can be written as 

dr 2 

ds 2 = -fir, 9)dt 2 + — — + E(r, 9)d9 2 + Kir, 9)d<\) 2 - 2H(r, 9)dtd<f>, (2.6) 
g{r,9) 

where /(r, 9),g(r, 9), S(r, 9), K(r, 9), H(r, 9) are defined below 

1 Q a 2 P sin 2 9 

/M) = n*^ - — 5— >. (2 - 7) 

9(r,6) = 9¥-, (2.8) 
P 

E(r,«) = ^, (2.9) 

The event horizons of this black hole can be calculated by putting 

— = 0, (2.12) 
9n 

which implies that 

g(r,6) = ^ = 0, (2.13) 
P 

or 

tt = 0,Q = 0, (2.14) 

which shows that in addition to the outer and inner horizons corresponding to the 
Kerr-Newman black hole horizons 

r ± = M± \/M 2 - a 2 , (2. 15) 
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we get two acceleration horizons [19, 201: r — - and r = — - 

° L ' ' a a cot 

We also define the function [9] 



Putting the values of f(r, 9), K(r, 9), H(r, 9) in this we get 



QPp 2 

(P(r 2 + a 2 ) 2 - Qa 2 sin 2 9)^ 



The angular velocity, for the above metric takes the form [9] 



Qh ~ K^-ry (2 - 18) 



Substituting H(r+,9), K(r+,6) we get 



a(P(r 2 + + a 2 )-Q(r + ) 
H Q(r + )a 2 sin 2 9 + P(r\ + a 2 ) 2 ' 1 ' ' 



Using the fact that Q(r + ) = 0, this can be written as 



n « = ^—-2- ( 2 - 2 °) 
r\ + a 2 

3. Tunneling of Dirac particles 

In order to study the tunneling of Dirac particles we solve the Dirac equation in the 
back ground of accelerating and rotating black holes. The covariant Dirac equation 
can be written as 

777 

i-f{Dp)* + -* = 0, (3.1) 

where we have 

D„ = d li + ft ft, (3.2) 

fy = X r ?b a ''A (3-3) 

and [7°, 7^] satisfy the commutation relations 

[7 a ,V , ] = -[7 /, ,7 a ] if [7 a ,7 /? ] = if a = /3. (3.4) 
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Here m is the mass of the fermions and 7^ matrices satisfy [7", 7^] = 2g^ u I, (I 
is the identity matrix). For fermion tunneling radiation, it is important to choose 
appropriate 7^ matrices. We take 



PQp 2 V p V P' 

*_ P^7 2 a(P(r 2 + a 2 ) - Q) 7 ° (g Q) 

sm6^P(r 2 + a 2 ) 2 -Qa 2 sin 2 9 ^/F{r, 6){P{r 2 + a 2 ) 2 - Qa 2 sin 2 9) ' 



Here 













<7 2 







(3.7) 



(3.8) 



The cr*(i = 1,2, 3) are the Pauli sigma matrices given by 















k-l 





(3.9) 



/I 

a 3 = -1 I . (3.10) 

In order to solve the Dirac equation we assume the following ansatz involving ar- 
bitrary functions of the coordinates corresponding to the spin-up and spin-down 
solutions of this equation 
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Here ^ and £4. are the eigenvectors of a 3 , and J t and I± denote the actions of the 
emitted spin- up and spin-down particles, respectively. On using Eq. (3.4) the Dirac 
equation takes the form 

771 

(t^dt + L^dr + L^de + i-y%)V + -jr* = 0. (3.13) 

Now, we substitute the above ansatz into the Dirac equation and compute it term by 
term. Dividing by the exponential term and neglecting the terms with h we obtain 
the following four equations 



= -B[ - 1 d t I t + 4 / ^-d r I t 



a(P(r 2 + a 2 )-Q) . „ , 

+ ^ — — 5— Vt + Am, 3.14 

V^FM)(P(r 2 + a 2 ) 2 - Qa 2 sin 2 0) rj ' 

= -B[J^fd e I t + ^ dM (3.15) 

V P sin e(y/P(r 2 + a 2 ) 2 - Qa 2 sin 2 6) * T 



= A[ ; 1 9t/ t - J^-dM 



V^FM)(P(r 2 + a 2 ) 2 - Qa 2 sin 2 0) ?J V ' 

= A [J^deh + = Vtl' ( 3 - 17 ) 

l V P 2 sin V^ 2 + a 2 ) 2 -Qa 2 sin 2 £ * f 

Taking into account the symmetries of the spacetime we assume the action to be of 
the form 

I t = -Et + J(f) + W(r,6). (3.18) 

Here E and J denote the energy and angular momentum of the emitted particle. We 
do calculations for the spin-up case; the spin-down case is similar. Inserting this into 
the above four equations, we get 



= -B[ - = + J^W'(r,e) 
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+ a(P(r 2 + a 2 )-g) ^ 

v/F(r, 0)(P(r 2 + a 2 ) 2 - Qa 2 sin 2 6) 

= B[J^fw e (r,e) + — J], (3.20) 

V P sm9{^P{r 2 + a 2 ) 2 -Qa 2 sin 2 

= X[-^=-./^W(r,«) 

wi V ^ 

a(P(r» + .')-<» j, + Bm< (3 . 21) 

v/F(r, fl)(P(r 2 + a 2 ) 2 - Qa 2 sin 2 9) 

= A[j¥fw e (r,0) + - J]. (3.22) 

V P sm9^/P{r 2 + a 2 ) 2 -Qa 2 sin 2 

Expanding Eq. (2.8) in Taylor's series and neglecting the higher powers we get 

g(r, 9) = g(r+, 9) + (r - r + )^ r (r + , 0). (3.23) 
Noting that at the horizon g(r + , 9) = 0, and its derivative at the horizon is 

2(l-qr + cosfl) 2 (r + -M)(l-q 2 r 2 ) 
^ = r 2 + a 2 cos 2 fl ' (3 ' 24) 

we obtain 

r+ + cr cos z 6* 

Similarly, noting that at the horizon F(r + ,9) = 0, from Eq. (2.17) we obtain 



F[r,9) = (r -r+)( 2 ; _^ 9r>9! ). (3.26) 



,2(r 2 + a 2 cos 2 0)(r+ - M)(l - a 2 r\ ) . 
(r 2 + a 2 ) 2 ^ 2 

Now expanding Eqs. (3.19) - (3.22) near the black hole horizon and using Eqs. (3.25) 
and (3.26) we get 



—E 



= -B[ = == + v / (r-r + ) 5r (r + ^)^'(r, 0) 

v/(r-r + )F r (r + ,0) 

I a(P(r i + a*)-Q(r + )) 

y/(r - r+)F r (r + , 9)(P(r 2 + + a 2 ) 2 - Q(r+)a 2 sin 2 0) 
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^ ^ r +'^ sm6^P(r 2 + + a 2 ) 2 - Q(r + )a 2 sin 2 



- ^ v / (r -; )fi ( r+ ,,) - V<r-r +)9r (r + , W r, S ), 

+ — J] + flm, (3.29) 

v/(r - r + )F r (r + , #)(P(r 2 + a 2 ) 2 - Q(r+)a 2 sin 2 9) 

= -J^-f.^ + «r + ,mr.,*) (3 , 0) 

V ^^ r +'^ sin9^P(r 2 + + a 2 ) 2 - Q(r + )a 2 sin 2 

We neglect the equations that depend upon "0", as their contribution to the total 
tunneling rate is canceled out, and retain only the radial equations (Eqs. (3.27) and 
(3.29)) which on using Eq. (2.19) take the form 



-E + Vl H J 



= -B\ ,J =f + V(r-r + )g r (r + ,9)W\r, 9)] + Am, (3.31) 

y/(r -r+)F r (r+,0) 

= A[ ~ E + QhJ == - y/(r-r + )g r (r + ,e)W'(r, 9)] + Bm. (3.32) 
y/(r -r + )F r (r + ,e) 

At the horizon we can further seperate W(r, 9) as 

W(r,6) = W(r) + &(9) (3.33) 
In the massless case there exist two possible solutions 

5 = 0, W'ir) = W'Jr) = (E _-p H J) , v 

+ y/(r-r + )F r (r + ,9)y/(r-r + )g r (r + ,ey 

and 

A = 0, W(r) = W'_{r) = (-£ +j Wj (3 35) 

V ( r - r +) ir r(r+, 0) V(r - r+)g r (r + , 9) 

Substituting the values of F r (r + ,9) and g r (r + ,9) the above equations become 

^ +(rj " 2(r-r + )(r+-M)(l-a 2 r 2 )' ^ 6b) 

, (-£ + lW)(r 2 +a 2 ) 

^- (r) ~ 2{r-r + ){r + -M){l-a 2 r 2 + Y (3 ' 37) 
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Here prime denotes the derivative with respect to r and +/— correspond to the 
outgoing/incoming solutions. For finding the value of W(r) we integrate the above 
result 

W (r) - f (E-n H J)(r*+a*) 

W+[r) ~ J 2(r - r + )(r + - M)(l - a*riy ^ 
Integrating around the pole r = r + this gives 

(E-n H J)(rl + a 2 ) 

Dropping the subscript we write 

W{T) = nL 2(r + -M)(l-a^ + y (3 ' 40) 

ImW = (3 . 41) 

2(r+ - M)(l - a 2 r\) V ; 

So the tunneling probabilities of fermions are 

Pemission = exp[-2/mi] = exp[-2(ImW+ + ImG)], (3.42) 

Pabsorption = exp[-2JraJ] = exp[-2(/miy_ + ImQ)]. (3.43) 

Since ImW + = —ImW_, we see that the total probability that the particle tunnels 
from inside the event horizon to outside is 

' Pem ™ = exp[-AImW+], (3.44) 



P 

1 absorption 

or 

r = g pff< g -"- J )(r Ur 2) l- (3-45) 
L l (r+ - M)(l - a 2 r\) J v 1 

Comparing this with Y = exp[— (3E] where j3 = 1/T/j we find that the Hawking 
temperature [6, 7] is given by 

Th = ^iL^ri) , (3 , 6) 

27r (r+ + a 2 ) 

where r + is given by Eq. (2.15). If we put acceleration equal to zero in formulae (3.45) 
and (3.46), they reduce to the tunneling probability and temperature of the Kerr 
black hole [8, 23]. Similarly, setting rotation equal to zero will recover expressions 
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for the Schwarzschild black hole. Comparing with the Kerr black hole, we note that 
the effect of acceleration is that it increases the temperature. 

From Eqs. (3.42)-(3.45) it appears that for some value of E and J the probabil- 
ities may become larger than 1 and hence violate unitarity. However, this is not the 
case because apart from the spatial contribution there is also a contribution to the 
imaginary part Im(EAt), from the temporal part of the action [24, 25]. This shifts 
the time by an imaginary amount which contributes both to P em ission and P a bsorption 
and results in a correct value of T. If we do not take this contribution into account, 
we will obtain the Hawking temperature twice as large as the actual value [26, 27]. 
Further, we note that Pabsorption will actually be equal to 1, because the trajectories 
of incoming particles do not face any barrier. This is, in fact, taken care of if the 
temporal contribution is also taken into account [24, 25], and we obtain the correct 
value of the tunneling probability. 

For the massive case Eqs. (3.27) and (3.29) do not decouple. We eliminate the 
function W'(r, 6) from these two equations by multiplying Eq. (3.29) by B and Eq. 
(3.27) by A and subtracting to yield 



A = -(E- m H ) ± y/(E - .m H f + m 2 F r (r + , 9)(r - 
P> m^J F r (r+, 9){r — r + ) 

In the limit r — > r + the two roots give either ^ — > or — > — oo, i.e. either A — > 
or B ->■ 0. For A ->■ we find the value of m from Eq. (3.32) 

m = ^( ~ E + ™^ = - V9r(r + ,e)(r-r + )W'(r)). (3.48) 
B y/F r (r+,6)(r - r+) 

Putting in Eq. (3.31) and simplifying we get 

Wr(r , o) = »-» = (E-m^ + xm 

+ V / ^(r + ,%,.(r + ,9)(r - r + )(l - A 2 /B 2 ) 

Integrating with respect to r as done before we finally get 



(E -Q H J)(r 2 + + a 2 ) 
v 1 2(r + - M)(l - a 2 r\y v 1 

(E-n H J)(rl + a 2 ) 

ImW = 71— — 7TTTT- 3.51 

2(r+-M)(l -a 2 r 2 + ) K ' 



For B — > we simply get 
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We note that we obtain the same tunneling probabilities as before, and hence the 
same temperature. This is because near the black hole horizon massive particles 
behave as massless. 



4. The acceleration horizon 

As mentioned earlier the black holes under consideration have an acceleration horizon 
at r a — 1/a also, apart from the rotation horizons. The functions F r (r, 6) and g r (r, 6) 
in this case will be 

_ (r 2 + a 2 cos 2 6)(a 2 -2Mr a + r 2 a )(-2r a a 2 ) 
r[ a ' ] (r 2 a + a 2 ) 2 il 2 (r a ,6) ' 

{l-cos9) 2 {a 2 a 2 -2Ma + l){-2a) 
9 ^ U >- (1 + aVcos 2 ^) ' 1 ] 

and from Eqs. (3.31) and (3.32) we see that the massless case gives rise to two 
possible solutions 

B = o, W \r) = W' + {r) = ^^f"^ w m - ^ 

yj(r - r a )F r {r a ,9)^/{r - r a )g r (r a ,9) 

and 

A = 0, W'(r) = W'_{r) = {-E + Q hJ} ^ 

y/(r -r a )F r (r a ,6)^(r -r a )g r (r a ,6) 

Putting values of the functions F r (r + ,9) and g r (r + ,9) in the above equations and 
integrating as before we obtain 

Similarly 

W ir)= -^E-n H J)(l + a 2 a 2 ) 
~ V } 2a(a 2 a 2 -2Ma + l) ' 1 ' ' 

Proceeding as before to find the tunneling probability, V = exp[— f3E], the resulting 
Hawking temperature at the acceleration horizon comes out to be 

a (a 2 a 2 - 2Ma + 1) 
H 27r(l + a 2 a 2 ) ' 1 ] 
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5. Calculation of the action 



Here we will work out the action explicitly. We have already seen that r- and 9- 
dependence decouples in Eqs. (3.27) - (3.30) near the black hole horizon. This 
allows us to separate the action as 

W(r,9) = R(r) + e(9). (5.1) 

Using this separation in Eqs. (3.31) and (3.32) we get 

= -B( - E + QhJ + ^(r-r + )g r (r + ,9)R'(r)) + Am, (5.2) 
y/(r - r + )F r (r+, 9) 

= A( - E + QhJ - ,/(r-r + )g r (r + ,e)R'(r)) + Bm. (5.3) 
y/(r - r + )F r (r+, 9) 

If m — 0, from Eq. (5.2) we get 

B = 0, R'(r) = R' + (r) = E j-Q. H J ^ 

V (r - r + )F r (r + , 0) v(r - r + )g r (r + , 9) 

Similarly from Eq. (5.3) we get 

A = 0, R'{r) = R' (r) = - £ + £l H J , v 

y/(r-r + )F r (r + ,6)y/(r-r + )g r (r + ,6) 

Substituting the values of F r (r + ,9) and g r (r + ,9) in these equations gives 



{E-n H J)(r 2 + + a 2 
2(r - r+)(r + - M)(l - a 2 r£) 

+ lW)(r 2 +a 2 ) 
2(r-r + )(r + -M)(l-a 2 r 2 ) 

Integrating Eqs. (5.6) and (5.7) we get 



K(r) = „,_ k tT"" J ^: V r (5-6) 
tM = ^TT^^^Lv (5-7) 



(E -tt H J)(ri + a 2 ) . . . . 

i2 (r) = -(^-^)K + ^ 2 ) ln(r _ r+) (5 9) 

[) 2(r + -M)(l-a 2 r 2 ) 1 +J ' 1 j 

For the massive case from Eq. (5.2) we get 

R> (r) = , Am + E - a » J (5.10) 

B^-^Sr (r + ,9) (r-r + ) y/g r (r + ,9) F r (r + 9) 
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where A and B are functions of (t, r,9,(p). After integrating it with respect to r and 
using the values of g r (r+, 9) and F r (r+, 0) we get for the outgoing particle 

/" Am (E - Q H J) (rl + a 2 ) , / x / 
R+ (r) = / =dr + ^ " ^ + - /. In r - r + . (5.11) 

Similarly from Eq. (5.3), for the incoming particle we get 

R- r = / ~ 777 , 9 \ In (r - r+) . 5.12 

W 7 V(r-r+)0r(r + ,0) 2 (r+ — M) (1 — a 2 r 2 ) V +; V > 

Now, we come to Eqs. (3.28) and (3.30) which take the form 



= -BU^-Se + = J), (5.13) 

'" P 2 smO{^P{r 2 + a 2 ) 2 -Qa 2 sin 2 9 

= -A(J^e e + = J). (5.14) 

1 P 2 sin 9{^P{r 2 + a 2 ) 2 -Qa 2 sin 2 #) 2 



Note that both the equations give 



0, = — ; Lp2j (5.15) 

sin 9VPy /p ( r2 + a 2 ) 2 - Qa 2 sin 2 9 



or at horizon 



(r+ + r)Psint) 
Substituting the values of p and P, and integrating we get 

Bm = -tJ f (rj + a 2 cos 2 9)d9 

U (r 2 +a 2 ) 7 sin0[l - 2aMcos# + a 2 a 2 cos 2 ' l ' J 

Using the method of partial fractions we evaluate the integral and finally get 



(9) = Li ln(l + cos 9) - L 2 ln(l - cos 9) + L 3 ln(l - 2aM cos 6> + a 2 a 2 cos 2 



aa 2 cos 9-M - y/W — a 2 \ , , / aa 2 cos 9-M - y/M 



+L 4 In - L 5 In . 

\aa 2 cos9 - M + y/M 2 - a 2 J \aa 2 cos9 — M + \JM 2 — a 2 



v-a 2 



where L { are 
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u 

2(1 + 2aM + a 2 a 2 Y 

U 

2(1 -2«M + aV)' 
iJaM 



(1 - 2aM + a 2 a 2 )(l + 2aM + a 2 a 2 ) ' 

iJa 2 [a 2 (1 + a 2 a 2 ) - 2M 2 ] 

2aVM 2 - a 2 (l - 2aM + a 2 a 2 )(l + 2aM + a 2 a 2 ) ' 
iJa 2 



2aVM 2 - a 2 (r 2 + a 2 ) 
This determines the action completely. 



6. Conclusion 

For studying Hawking radiations from black holes different types of approaches have 
been adopted in the literature. This has been done by using the Newman- Penrose 
formalism and the so-called Hamilton- Jacobi method. In particular, tunneling of 
Dirac particles has been studied for the Kerr and Schwarzschild black holes. We 
have extended this semi-classical approach to study a large class of black holes which 
include those with acceleration and rotation. We obtained the tunneling probability 
for the incoming and outgoing particles and correctly recover the Hawking tempera- 
ture. We have explicitly calculated the action as well. An indication of the generality 
of our results is that, in appropriate limits, they reduce to those for the Kerr and 
Schwarzschild black holes. 

Acknowledgments 

We are thankful to Douglas Singleton for some useful comments. 
References 

[1] S. W. Hawking, Commun. Math. Phys. 43 (1975) 199. 
[2] S. W. Hawking, Nature 248 (1989) 30. 

[3] M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85 (2000) 5042 [hep-th/9907001]. 



-14- 



[4] P. Kraus and F. Wilczek, Nucl. Phys. B 433 (1995) 403. 

[5] M. K. Parikh, Gen. Rel. Grav. 36 (2004) 2419 [hep-th/0405160]. 

[6] K. Srinivasan and T. Padmanabhan, Phys. Rev. D 60 (1999) 24007. 

[7] S. Shankaranarayanan, T. Padmanabhan and K. Srinivasan, Class. Quant. Grav. 
19 (2002) 2671. 

[8] R. Kerner and R. B. Mann, Phys. Rev. D 73 (2006) 104010. 

[9] R. Kerner and R. B. Mann, Class. Quant. Grav. 25 (2008) 095014. 

[10] R. Kerner and R. B. Mann, Phys. Lett. B 665 (2008) 277. 

[11] R. Li, J. R. Ren and S. W. Wei, Class. Quant. Grav. 25 (2008) 125016. 

[12] D. Y. Chen, Q. Q. Jiang and X. T. Zu, Class. Quant. Grav. 25 (2008) 205022. 

[13] X. X. Zeng and S. Z. Yang, Gen. Rel. Grav. 40 (2008) 2107. 

[14] C. Ding and J. Jing, Class. Quant. Grav. 27 (2010) 035004. 

[15] J. Yang and S. Z. Yang, J. Geom. Phys. 60 (2010) 986. 

[16] R. Li and J. R. Ren, Phys. Lett. B 661 (2008) 370. 

[17] R. Li, S. Li and J. R. Ren, Class. Quantum Grav. 27 (2010) 155011. 

[18] J. Podolsky and H. Kadlecova, Class. Quantum Grav. 26 (2009) 105007. 

[19] J. B. Griffiths and J. Podolsky, Class. Quantum Grav. 22 (2005) 3467. 

[20] J. Podolsky and J. B. Griffiths, Phys. Rev. D 73 (2006) 044018. 

[21] O. J. C. Dias and J. P. S. Lemos, Phys. Rev. D 67 (2003) 064001. 

[22] O. J. C. Dias and J. P. S. Lemos, Phys. Rev. D 67 (2003) 084018. 

[23] M. Bilal and K. Saifullah, Thermodynamics of accelerating and rotating black holes 
[arXiv: 1010.5575]. 

[24] V. Akhmedova, T. Pilling, A. de Gill and D. Singleton, Phys. Lett. B 666 (2008) 
269. 

[25] V. Akhmedova, T. Pilling and D. Singleton, Int. J. Mod. Phys. D 17 (2008) 2453. 

[26] E. T. Akhmedov, V. Akhmedova and D. Singleton, Phys. Lett. B 642 (2006) 124. 

[27] T. Pilling, Phys. Lett. B 660 (2008) 402. 



- 15 - 



